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The creation of collective many-body quantum states from a two-dimensional lattice gas of atoms 
is studied. Our approach relies on the van-der-Waals interaction that is present between alkali metal 
atoms when laser excited to high-lying Rydberg s-states. We focus on a regime in which the laser 
driving is strong compared to the interaction between Rydberg atoms. Here energetically low-lying 
many-particle states can be calculated approximately from a quadratic Hamiltonian. The potential 
usefulness of these states as a resource for the creation of deterministic single-photon sources is 
illustrated. The properties of these photon states are determined from the interplay between the 
particular geometry of the lattice and the interatomic spacing. 
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I. INTRODUCTION 



The quantum interface between light and an atomic 
ensemble has attracted much attention during the past 
decade Jl:. The reason is rooted in the wide range of 
possible applications that rely on the coherent coupling 
between the two systems: Implementation of quantum 
information processing protocols [2j [3], slow-light and 
electromagnetically induced transparency [3HE] as well 
as the creation of deterministic and manipulable photon 
sources [9HT2]. In order to create these photon sources, 
one must be able to create entangled atomic many-body 
states and map them efficiently into the desired photonic 
states [T5HTB] . Ultracold atoms represent the ideal tool 
for creation these atomic states under very clean and 
well-defined conditions. This is due to the advanced tech- 
niques that are nowadays available for the trapping of ul- 
tracold atoms and for tailoring the interactions between 
them [17] ■ Highly excited (so-called Rydberg) states [18] 
have proven to be especially useful in this context due 
to the strong state-dependent interaction between them, 
allowing the entanglement of atoms separated by several 
micrometers [19H23) . 

In a recent work a laser-driven gas of Rydberg atoms 
trapped on a one-dimensional ring lattice has been 
proposed as a resource to create collective excitations 
[23] [53] • It was shown that these excitations could be 
converted into few photons following a mapping scheme 
proposed in Re fs. [T3R15] [26] . As a result non-classical 
states of light such as single-photon sources and entan- 
gled pairs of photons could be created [16] whose prop- 
erties where determined by the interplay of the ring ge- 
ometry and the nature of the delocalized atomic excita- 
tions. In this work, we extend this study to a gas of 



atoms trapped in two-dimensional medium-scale square 
and triangular lattices. We study thoroughly the prop- 
erties of collective many-body states in these two lattice 
types and outline methods of selective excitation. We 
find that one can create single-photon sources in which 
the photon is emitted in a superposition of a few well- 
focused beams. The direction of these beams is mainly 
determined by the geometry and spacing of the lattice. 

The paper is structured as follows. In Section [IT] we 
introduce the system we have in mind and present a 
derivation of the Hamiltonian that governs its dynam- 
ics. In Section |III| we discuss the diagonalization of this 
Hamiltonian and explain the approximations considered 
to do so. The symmetry properties of the many-particle 
eigenstates are analyzed in Section IV together with a 



proposal of how to excite experimentally these states. 
We explain how to map these atomic collective excita- 
tions into non-classical states of light and discuss their 
properties in Section [V| Finally, in Section VI we discuss 
the validity of the approximations made and take into 
account disorder caused by quantum uncertainty and/or 
temperature. We conclude with a summary and outlook 
in Section [VHl 
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II. THE SYSTEM AND ITS HAMILTONIAN 

In the setup we have in mind a gas of atoms is trapped 
in a regular two-dimensional lattice with N sites and a 
lattice spacing a of the order of a few micrometers [27l - 
f5Uj (see Fig. [lji) . We consider the limit of a deep lattice 
where the vibrational states within each site are well- 
approximated by the eigenstates of a harmonic oscilla- 
tor and tunneling between the sites is absent. The lat- 
tice shall be uniformly occupied with one atom per site 
(Mott-insulator state) and we assume zero temperature 
such that each atom populates the ground state of its po- 
tential well whose spatial width a is much smaller than 
the interparticle separation a. 
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The internal (electronic) degree of freedom of each 
atom is approximated by a two-level system. The ground 
state \g) is coupled to a Rydberg ns-state |r) by means of 
a laser with Rabi frequency £1 and detuning A as shown 
in Fig. [T]d. For the moment, we do not focus on the two 
auxiliary levels \s) and \a) also shown in this figure. They 
will be used for the excitation of the many-body states 
and the mapping of these states into light, that will be 
explained thoroughly in Sections IV and|Vj respectively. 
Within the two-level approximation we can identify each 
atom as a spin- 1/2 degree of freedom where \g) k = \\) k 
and \r) k = |t) fc [?TH55] . We can thus employ the Pauli 



matrices (ai k \ a y k ' and of"') to formulate the Hamil- 
tonian governing the excitation dynamics of our system. 
Using this the Hamiltonian describing the laser excita- 
tion process reads in the rotating wave approximation 



H h = 



N 

E 

fe=i 



[flajf 



(k) 



Ani- 



is the Rydberg number oper- 



where n k = 1/2 (l + a 

ator of the fc-th site. 

When the atoms are in the Rydberg state, they interact 
strongly via van-der-Waals interaction K,dw( r ) = C-6 x 
|r|~ , where r is the separation between the atoms and Cq 
is the van-der-Waals coefficient [351 ISZ] ■ The interaction 
Hamiltonian thus becomes 



H, 



V km n k n r , 



with the coefficients 



Vkm — 



r fe 



where r k denotes the position of the atom trapped in 
the fc-th site. Adding the interaction with the laser, the 
complete Hamiltonian of the system reads 



N 



k=l 



flai k) + An k 



E ^ 

m^k 



,n k n r 
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We assume the external dynamics of the atoms to be 
frozen. This assumption is well justified because the in- 
ternal dynamics take place on a much shorter time scale 
(typically of the order of a hundred nanoseconds versus 
milliseconds). We will also assume for the moment that 
the atoms are infinitely localized at the center of the 
traps, i.e. a/a — > 0. This very idealized situation can 
only be achieved approximately in experiment. We will 
investigate the effect of a finite width of the wave packets 
in Section fVll 



III. DIAGONALIZATION 

In this section, we will show how to obtain the eigen- 
values and eigenstates of the Hamiltonian and discuss 
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FIG. 1: a: Two-dimensional lattice in which the atoms are 
trapped. The width of the external wave function of each 
atom is given by a which is considered to be much smaller 
than the lattice spacing a. The largest interaction energy 
in the system is the one between nearest neighbors V nn . b: 
Internal level structure of each of the N atoms on the lattice. 
The two ground states \g) and \s) are coupled to the Rydberg 
state \r) by means of two classical lasers with Rabi frequencies 
Q and £72 , respectively. A third laser £7i drives the transition 
\g) — > \s). Finally, \g), \s) and an auxiliary state \a) form a 
lambda scheme where \s) is coupled off-resonantly to \a) and 
photons are emitted in the transition |o) — > \g). 



the approximations used to obtain them. Throughout 
this work, we will focus on the regime where the laser 
driving ft is the largest energy scale of the system. In ad- 
dition, the detuning of the laser |A| is taken to be much 
smaller than both the Rabi frequency (laser driving) and 
the interaction strength such that f2 » V nn ^> |A|. 



Holstein-Primakoff transformation 



The Hamiltonian ([I]) contains spin operators that obey 
anti-commutation and commutation relations when they 
belong to the same and different sites, respectively. 
Therefore, the underlying algebra is neither bosonic nor 
fermionic. In order to be able to solve the system at 
least approximately the spin operators will be expressed 
in terms of operators that obey a purely bosonic algebra. 
This is achieved by the Holstein-Primakoff transforma- 
tion [38] . 

We will make use of a slightly unconventional formu- 
lation of this transformation. As fl ^> V nn 3> |A|, the 
first term of Hamiltonian ([T]) is the dominant one and it 
is thus convenient to make it diagonal. For this reason, 
instead of the usual form of the Holstein-Primakoff trans- 
formation, we will use a rotated version of it such that 



„(*0 



a^\ This is achieved by means of the unitary 



operator U k = exp 



r -a {k) 
i y 



after which the Holstein- 



3 



Primakoff transformation becomes 



B. Constrained dynamics 



r (fc) 



= V 1 ~ a k a k a-k 



(k) 



(*n - „t 



a\Jl-a\a k (2) 



<4 fe) = 2a\a k - 1, 



where (aj.) are bosonic operators that create (annihi- 
late) non-interacting bosonic excitations within the sys- 
tem at the fc-th site. Considering that the eigenstates of 
the single-atom operator are 

1 



fcj ' 



the equivalent bosonic states are 



|0), and |+), 



U) fe = 4. 



ith zero and one bosonic excitation on 



site k, respectively. 

By making the substitution ([2]), the Hamiltonian ([!]) 
takes on a rather complicated form that contains square 
roots of operators, which are difficult to treat in practical 
calculations. In order to overcome this difficulty, we make 
an approximation: our study will be focused on many- 
body states that carry only a few excitations, i.e. whose 

total number of bosonic excitations Ab = 53& ^ a l°fc^ is 
much smaller than the number of sites N. In the sub- 
space spanned by these states - denoted by (...) - one has 
(a k a k } -C 1 for all k. Dropping the notation (...) in the 
following, one can make a Taylor expansion of the square 
root operator such that, after the normal ordering of the 
operators, one obtains 



1 — a k a k ~ 1 — a k a k 



We now substitute the operators of eqs. ([2| into Hamilto- 
nian Q and consistently with the condition (a£ a fc) <C 1 
keep only the first two terms of the expansion of the 
square root. Furthermore, we neglect all terms of the 
Hamiltonian that are of higher order than quadratic in 
a k . The Hamiltonian then assumes a quadratic form 



JV 



fc=i 

(4 + - 2 51 Vkm (4 + a > 

k—1 k^m 

7 5Z Vkm (4 a m + afeflm + 4 a ™ + a « 



(3) 



A 
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where the energy offset Eq is given by 



(4) 



A closer inspection of the Hamiltonian pi reveals that 
it is composed of three distinct parts {H = Eq + Hq + 
Hi + H2) which effectuate qualitatively different cou- 
plings between states with a fixed number of bosonic ex- 
citations. The first part is formed by those terms which, 
applied to a given state, do not change its number of 
excitations, i.e. 



N 

H = 2n^2 a{a k + V km [a\a m + a 



(5) 



k=l 



as they contain the same number of creation and anni- 
hilation operators. Second, there are a number of terms 
that create or annihilate a single bosonic excitation, 



A 

Fl =-2 



N 1 — , 

(4 + ak ) - 2 H Vkm ( c 



a h . 
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Finally, there is a third class of terms, which change the 
total number of bosonic excitations in the system by two: 



H 2 = \ Vkm (44* 



a k a r , 



(7) 
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FIG. 2: Scheme of the low-energy sector of the spectrum, 
formed by quasi-degenerate manifolds of states with the same 
number JVb of bosonic excitations separated by an energy gap 
of approximately 2Q. The arrows represent schematically the 
transitions driven by the different parts of the Hamiltonian 

Due to the strong laser driving (fi ^> V nn |A|), the 
first term of the Hamiltonian Hq which is proportional 
to f2 is clearly dominant and thus determines the coarse 
structure of the spectrum of the system. This coarse 
structure is formed by manifolds of quasi-degenerate 
states separated by energy gaps of approximately 2Q. 
As it is shown schematically in Fig. [2] each manifold 
is formed by states that contain the same total number 
Ab of bosonic excitations. Taking this into account, let 
us now analyze the rest of the terms of the Hamiltonian 
one by one. The second term in Hq does not change the 
total number of excitations of the state. Its action pro- 
duces couplings between states inside the same energy 
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manifold, whose strength is proportional to V nn . The 
Hamiltonians H\ and Hi create or destroy one and two 
excitations every time, respectively. Since the states in- 
volved in these transitions are separated energetically by 
approximately 2f2 and 40, respectively, they are strongly 
suppressed. In particular, the corresponding transition 
rates can be estimated by second order perturbation the- 
ory to be proportional to A 2 /£l and V nn2 /fl, which in 
our parameter regime (fi 3> |A|,y nn ) are very small 
quantities. As a consequence it can be affirmed that, to 
a very good degree of approximation, the Hamiltonian 
([5]) that drives the intra-manifold dynamics also governs 
the dynamics of the entire system. The validity of this 
approximation will be assessed in Section [VI] 



C. Eigenexcitations 



The manifold with Nb = 2 is formed by all combinations 
of two of the N eigenexcitations 



b\b\ |0) , 



So 



E 



(10) 



with j > i = 1 . . . N, so that there exist a total num- 
ber N(N + l)/2 of these doubly excited states. Follow- 
ing the same procedure, one could think of constructing 
states with higher numbers of excitations by applying 
more creation operators to the vacuum. However, one 
has to keep in mind that the Hamiltonian ([5]) gives a 
valid description of the system only if (a£afc) <C 1, i.e. in 
the low-occupation sector of the spectrum where only a 
few bosonic excitations compared to the total number of 
atoms are present. 
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We will now diagonalize the Hamiltonian Ho [eq. ([5 
Since it is quadratic this is achieved by means of t 
diagonalization of the matrix V, that essentially contains 
the interaction energies between each pair of atoms. This 
is achieved by a unitary transformation which reads 



J2 MktDiSijM}, 



where M k i represents the fc-th entry of the i-th eigenvec- 
tor of V and Di is its corresponding eigenvalue. Sub- 
stituting the last expression into ([5]), the diagonalized 
Hamiltonian becomes 



JV 



H = J2 e Mbi, 



i=l 



with the operators 



b = 



JV 

E 

fc=i 



M ki a\ 



and the eigenvalues e, = 20 -I- ^ for i = 1, . . . , N. Hence, 
the collective excitations of the system are formed by a 
superposition of all the bosonic excitations at, with the 
coefficients of the superposition given by the eigenvectors 
of V. 

We can now write down the eigenstates and corre- 
sponding energies of the system. The ground state is 



JV 



io> = n i°) 



(8) 



fe=i 



and its energy is E , given by eq. Q. The first excited 
manifold is spanned by the N states that possess one 
collective bosonic excitation, i.e. 



U) = b] |0) , E u =E a + e t 



(9) 



IV. MANY-PARTICLE STATES: SELECTIVE 
EXCITATION AND PROPERTIES 

Our aim is to selectively excite the many-particle states 
described in the previous section starting from the ex- 
perimentally relevant initial state in which all atoms are 
prepared in the electronic ground state 



init = 



(ii) 



To this purpose, we use the excitation scheme proposed 
in Refs. [53J[2S], that consists of two steps: 

First, we want to achieve the ground state of the 
Hamiltonian H$, |0) = IlfcLi |0)fc' wmcn contains zero 
bosonic excitations. To do so, we express the single-atom 
bosonic states in terms of the atomic ones, so that the 



ground state yields |0) = Ilfc=i \~)k- W e also introduce 
an additional internal single atom stable state \s) k (see 
Fig. [T]d). In practice \s) k could correspond to one of the 
states of the hyperfine ground state manifold which is 
different from \g) k . Using this state we can prepare the 
desired state |0) by means of the two subsequent laser 
pulses which are depicted in Fig. [TJd. We choose the first 
pulse, with Rabi frequency Qi, to be resonant with the 
single-atom transition \g) k — > |s) fe . After applying a 7r/2- 
pulse during a time n = 7r/(2£li), this laser is turned off 
and another one with Rabi frequency f2 2 and resonant 
with \s) k — > \r) k is shone on the system during a time 
T2 = 7r /VL2. This amounts to the sequence 



fc +*i s )j 



(12) 



-4 



V2 



=m->fc = i°> 



which results in the desired preparation of |0). Here we 
would like to remark that the strength of the second laser, 
given by its Rabi frequency has to be large enough 
to overcome the interaction energy between atoms in 
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the Rydberg state. The required strength will depend 
strongly on the particular geometry of the system. For 
example, in a square lattice each atom has four nearest 
neighbors, so that in this case the condition to be ac- 
complished is roughly O2 3> 4V nn , while in a triangular 
lattice it would be 3> 6^"" as there are six nearest 
neighbors for each site. 

Once the ground state of the Hamiltonian §5§ is pre- 
pared, the aim is to excite the previously described 
many-particle states (|9| and ( 10 ) that contain one and 
two bosonic excitations, respectively, thus belonging to 
higher energy manifolds. The idea proposed in Refs. 
[231 [35] is to induce transitions between states with 
AiVb = ±1 making use of the term of the Hamiltonian 
proportional to the detuning, i.e. the first term of 



k=l 



ak 



(13) 



which creates or annihilates one bosonic excitation at 
a time. The problem is that, as we discussed before, 
these transitions are highly suppressed due to the large 
value of the energy gap - approximately 20, see Fig. [2] 
- between the manifolds. To overcome this difficulty one 
can introduce a time-dependent detuning of the form 

A(t) = A cos (u>At), 

which will therefore act as a coherent oscillating source 
that couples the two many-body states considered when 
the oscillation frequency uj& matches the corresponding 
energy gap between them. The intensity of the transition 
between the two states considered is calculated in each 
case by working out the square of the transition matrix 
elements / = |/)| 2 , where \i) and |/) represent 

the initial and final state, respectively. 

Let us first consider the excitation from the ground 
state |0) to the first excited manifold of states (|9| that 
contain a single collective bosonic excitation. The inten- 
sity yields in this case 



j x (i) = Ko|fl- A |ii>| 2 = ^ 



N 



k=l 



(14) 



where the Rotating Wave Approximation (RWA) has 



been made and ti 



2VL 



is the corresponding en- 



ergy gap between the state and the ground state. 
Note that the RWA implies that the value of the energy 
gap - which is roughly proportional to Q - is much larger 
than the detuning, i.e. Ao <C f2. 

Following the same procedure, one can calculate the 
intensity of the transition between the one- and two- 
excitation states, that yield 

I 2 (i;jk) = |(1^ A |2^)| 2 (15) 

1 + 5 3k 



Note that this expression already tells us that the inten- 
sity of the transition between the two states will be zero 
unless one of the excitations in the two-boson state is 
the same as the single one originally present, i.e. unless 
i = j or i = k. Other selection rules that stem from the 
particular geometry of the system will be analyzed in the 
following section. 



A. Selection rules 

When calculating the intensity of the transitions be- 
tween the different manifolds following Eqs. (14) and 



(15), we observe that only for a few states they have a 
non-zero value (see Figs. 4] and [5]). The reason for this 
can be found considering the particular symmetries of the 
system. Let us consider here two examples of geometries: 
a square and an equilateral triangular lattice. 
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FIG. 3: Scheme of two lattices and their symmetries, a: The 
group dv that contains the symmetries of a square is formed 
by 8 elements: Rotations of the lattice by 0, n/2, w and 3n/2; 
flip about the x and y axes; flip about the diagonals u and v. 
b: The group Cz v contains the symmetries of an equilateral 
triangle, formed by 6 elements: Rotations of 0, 27r/3 and 47r/3 
and flip about the axes a, b and c that go through the vertices 
of the triangle perpendicular to the opposite side. 

The symmetry groups corresponding to the square and 
the equilateral triangle are C^ v and C^ v , respectively. 
These groups are formed by the following transforma- 
tions (see Fig. [3^i and b): 

a. Square 

• C*, rotations of 27rn/4, with n — 0, 1, 2, 3 

• F x ,F y , flip about the vertical and horizontal axes 

• F U ,F V , flip about the two main diagonals 

b. Equilateral triangle 

• C^, rotations of 27rn/3, with n = 0, 1, 2 

• F a , Ft,, F c , flip about the three axes trough the ver- 
tices 

One can observe that both the complete Hamiltonian H 
([lj and i?A ( 13 1 conserve all the symmetries correspond- 
ing to the system's geometry. As a consequence, if the 
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initial state is an eigenstate with respect to the previous 
symmetry operators, the time evolution will be restricted 
to the subspace spanned by the states with the same 
quantum number with respect to the cited operators. Ac- 
tually, both the experimentally relevant initial state |init) 
and the ground state |0) where no bosonic excitation is 
present, are eigenstates of the operators that represent 
all the rotations and flip operations that form the cor- 
responding groups with eigenvalue +1. This is true for 
both the square and the triangular lattice. The states 
that have eigenvalue +1 with respect to all the operators 
representing the group elements belong to the subspace 
A±, called totally symmetric. Thus, only the states that 
belong to this totally symmetric subspace can be accessed 
in the course of the time evolution under Hamiltonians 
H and Ha- 
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FIG. 4: Normalized absorption profile for the transition |0) — > 
Only the states that belong to the totally symmetric 
subspace Ai have a non-zero probability of being excited. The 
three states associated to the highest intensities are shown a: 
in a 7 x 7 square lattice and b: in a equilateral triangular 
lattice with 7 sites per side. The energies are given in units 
ofV nn . 

Keeping these selection rules in mind, let us now con- 
sider the excitation of the states with a single collective 
bosonic excitation from the ground state |0) and calcu- 
late the intensity of these transitions. An example is 
given in Fig. [4^, and b for a square lattice of N = 49 
sites and a triangular one with N = 28, respectively. 
Here we observe that, as a consequence of the selection 
rules, the intensity is always zero except for a very re- 
stricted set of final states. These states belong to the 
subspace A±, i.e. they are symmetric (eigenvalue +1) 
with respect to all the symmetry operations correspond- 
ing to the particular geometry of the lattice. Finally, we 
also observe that the state with highest eigenenergy is 



the one with the largest intensity. This is observed for 
all sizes of the system (determined by N) that have been 
studied and for both geometries considered. This state 
reads = J^k ^kN% |0), with M^n > 0, and is the 
closest state to a uniformly shared collective excitation, 
i.e. a spin wave without spatially varying phase factors. 
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|liv) 



5: Normalized absorption profile for the transition 
-> |2jj). Only the states that belong to the totally sym- 



metric subspace A\ have a non-zero probability of being ex- 
cited. The density of bosons of the two states associated to 
the highest intensities are shown a: in a 7 x 7 square lattice 
and b: in an equilateral triangular lattice with 7 sites per 
side. The energies are given in units of V nn . 



Let us now discuss briefly the excitation of the states 
that carry two bosonic excitations. In particular, let us 
assume that we have already excited the state 1 1 pj) which 
corresponds to the largest intensity h(N) and calcu- 
late the transitions from that state to the doubly-excited 
manifold. As we already obtained in Eq. ( 15 ), the inten- 



sity is zero unless one of the two excitations corresponds 
to the initial one, e.g. in this case fejy. This corresponds 
to only N possible states out of the total N(N + l)/2 
available. Moreover, the selection rules will reduce the 
number of accessible states even more due to the fact 
that only the states that belong to A\ (totally symmet- 
ric states) can be accessed from |1jv), which also belongs 
to A\. In particular, a state that carries two excitations 
only belongs to the symmetry subset A\ when the two ex- 
citations that form it belong to the same symmetry sub- 
space. As a consequence, the number of doubly-excited 
states with non-zero transition intensity is the same as 
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the number of singly excited ones with Ji ^ 0, which 
can also be directly inferred from Eq. (15). To under- 
stand this better, let us study an example. For a 3 x 3 
[N = 9) square lattice there are only three eigenexcita- 
tions that belong to A\ , which are b\ , b\ and b\ . Let us 
consider that we are interested in knowing which states 
can be then accessed from the state |lg) = b\ |0). In prin- 
ciple any combination of the eigenexcitations bjb k with 
j > k = 1 . . . 9 should be accessible. However, we know 
from the result (15) that only if j =9 or k = 9 the 



intensity of the transition is non-zero. Finally, the selec- 
tion rules impose that only the combinations of the three 
eigenexcitations that belong to Ax, i.e. |2ig), |2sg) and 
|2gg) can be accessed from |lg) with I2 ^ 0. In Fig. [5^, 
and b we show the corresponding intensities for the same 
two systems considered before, i.e. 7x7 square lattice 
and equilateral triangle with 7 sites on each side, respec- 
tively. In this case, in the insets of the figures show the 
bosonic density in each site n, i.e. {2j k \a\ l a n \2j k ) , for 
the two states corresponding to the highest intensity. 



B. Mapping to a stable state 

Once the many-particle atomic states have been ex- 
cited, we need to map them into a stable configuration. 
This is necessary due to the limited lifetime of the Ry- 
dberg states and if one wants to manipulate the corre- 
sponding state further. To this end we make use of the 
stable hyperfine states introduced in the previous section, 
\s) k (see Fig. [ljb) . Let us remind that each of the single- 
atom bosonic states |l) fc , |0) fe is equivalent to the atomic 
single-atom one |+), and |— ). , respectively (see Section 



III A). As a consequence, by inverting the sequence of 
laser pulses given in (|12[), we achieve the mapping 



|0> fc ^lff) ? 



ID 



so that the excitations are stored now in the two hy- 
perfine ground states. Expressed in terms of the atomic 
states, the many-body states studied and given by ([9| 
and (10) become 



(16) 



k=l 



and 



N 



*«) = E M kl M l3 a^^ g |init) , (17) 



fcj=i 



respectively, with a 



(k) 



(si- 



V. APPLICATION: SINGLE PHOTON 
SOURCES 

So far, we found the low-lying eigenexcitations in the 
lattice and described the properties of the many-body 



states that can be accessed experimentally by controlling 
the laser parameters. In this section, we will show how 
to map the collective atomic excitations into quantum 
states of light [T3HT5] . 



A. General aspects of the atom-photon mapping 

The mapping from the atomic excitations encoded in 
the two hyperfine ground states \g) k and \s) k to non- 
classical states of light has been thoroughly investigated 
in [T3J [T5]. In the following we will briefly outline the 
basic idea. 

We consider the level scheme shown in the lower part 
of Fig. [Tjb where the two ground states \g) k and \s) k to- 
gether with an auxiliary state \a) k form a lambda scheme. 
The storage state \s) k is coupled off-resonantly to the 
state \a) k via a classical laser field with Rabi frequency 
and detuning and Al respectively. Eventually, a pho- 
ton is emitted to the radiation field by means of the tran- 
sition \a) k — > \g) k - Note that we do not consider the de- 
cay from the state \a) k back to \s) k , a condition that can 
be accomplished by an appropriate choice of the atomic 
levels [13]. 

This setup allows us to establish a direct mapping be- 
tween any given atomic state - encoded in \g) k and \s) k 
- and its corresponding photonic one after a time much 
longer than the lifetime r = T^ 1 of the intermediate state 



\a) k with r being the corresponding decay rate. In par- 
ticular, one can show that a single atomic excitation is 
mapped to a single photon state as 

where a qi/ creates a photon with momentum q and po- 
larization v. The details of the coupling are contained in 
the function <?fc qi , (t) , that - in the limit of t 3> r - is given 
by 



N 



71,7=1 



-1 

n/v 

ILU — K 7i 



. (18) 

Let us clarify the meaning of each term appearing in the 
previous expression. First, kL is the momentum of the 
laser and lu = |q| /c is the energy of the emitted photon, 
which also appears in the coefficient 



^ ~ Al V 2e V ga qi " 

Here, V is the quantization volume, eo the vacuum per- 
mitivity, d ga the dipole operator of the \g) k — > \a) k tran- 
sition and e q „ the polarization vector (q • e qt , = 0). Fi- 



nally, in eq. ( 18 ) Xn and K n represent respectively the 



eigenvectors and eigenvalues of the operator that gov- 
erns the atomic dynamics in the process of light emission 
from an ensemble of atoms. This operator accounts for 



light scattering at multiple atoms as well as interatomic 
dipole-dipole interactions and depends mainly on the rel- 
ative orientation of the atomic transition dipole moments 
and the ratio between the interparticle separation and the 
wavelength of the laser, o/Al [2B"] - 



B. Creation of single photons 

We will now use this mapping to obtain the single- 
photon states that are created when using the atomic 
states ( 16 1 that contain one collective excitation as a re- 



source. The produced single-photon states are 



JV 



®i ) = E E M ki9k^(t)a{ v | vac) , 

fc=l qy 

where |vac) represents the photon vacuum. In order to 
characterize the properties of the photonic states, we will 
study the angular photonic distribution, i.e. 

V f°° 



where 



is the number of photons with mo- 



mentum q and polarization v. This quantity provides 
us with the average photon number in the direction of 
emission parameterized by 9 and (f>. 

Let us now specify the parameters of our setup. We 
consider the lattice to be lying in the xy-pl&ne. The 
dipoles of the transitions \g) k — > \a) k are all aligned and 
oriented perpendicular to the lattice, i.e. d gQ |j z. This 
can be ensured by an appropriate choice of the atomic 
levels. Finally, we consider the situation where the inci- 
dent direction of the laser momentum Icl is also perpen- 
dicular to the lattice plane, i.e. || z. The only degree 
of freedom left is - apart from the specific geometry of 
the lattice - the ratio between the interparticle separa- 
tion given by the lattice spacing a and the wavelength of 
the laser Al ■ We will see in the following how the change 
of this ratio modifies dramatically the properties of the 
emitted photon. 

In Figs. [6] and [7] we show two examples of these single 
photon sources. In both cases we have chosen to map to a 
photonic state the atomic many-particle state with a sin- 



gle collective excitation with the highest energy 

which also correspond to the largest transition intensities 
for each geometry (see Figs. [4^, and b) . The correspond- 
ing states are represented in Figs. [6^, and [7^ for a square 
and a triangular lattice, respectively. The red arrow rep- 
resents the direction of the incident mapping laser. For 
each geometry we show the angular distribution of the 
emitted photon for two different values of the ratio o/Al- 
In the case of a square lattice and for o/Al = 0.9 (see 
Fig. [6Jd) we see that the photon is emitted in a super- 
position of four directions, each one perpendicular to the 
sides of the square lattice. For the same ratio and a 






FIG. 6: a: The atomic state * 



(i) 



in the case of a 7 x 7 



lattice is mapped into a single photon state by means of a 
laser whose direction is given by the red arrow, b: When 
(i/Al = 0.9, the photon is emitted in a superposition of four 
beams that are perpendicular to the sides of the square lattice, 
c: The angular distribution is strikingly different in the case 
os/Al = 0.25, where the emission occurs into two pyramidal- 
shaped beams with square base, as it can be observed in d. 
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FIG. 7: a: The atomic state |\t^ ) in the case of a triangu- 
lar lattice of side 7 is mapped into a single photon state by 
means of a laser whose direction is given by the red arrow, 
b: When a/X^ — 0.9, the photon is emitted in a superposi- 
tion of three beams that are perpendicular to the sides of the 
triangular lattice, c: The angular distribution is strikingly 
different in the case a/X^ — 0.4, where the emission occurs 
into two pyramidal-shaped beams with triangular base, as it 
can be observed in d. 



triangular lattice (Fig. [7Jd) , the photon is emitted in a 
superposition of three directions, again perpendicular to 
the sides of the lattice. With this example we see clearly 
how the geometry of the lattice plays an important role 
in the properties of the photonic state. These states are 
particularly interesting since the photon in this case is 
not only emitted in a highly directional form but also 
perpendicularly with respect to the incident laser, which 
would avoid noise problems in the photon detection pro- 
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In order to illustrate how the ratio 0,/Al also affects the 
properties of the state, we show in Fig. [6J: the mapping 
of the same atomic state in the square lattice but this 
time with <j/Al = 0.25. We observe that in this case 
the photon is emitted in a superposition of directions 
that form two cone-shaped beams that still conserve the 
information of the geometry of the lattice as can be seen 
in Fig. [6]i. The same numerical experiment has been 
done in the case of the triangular lattice in Figs. [7J: and 
d for o/Al = 0.4 with a similar result. 



VI. VALIDITY OF THE APPROXIMATIONS 
AND ROBUSTNESS OF THE SOLUTIONS 

So far, we have considered a very idealized situation. In 
particular, we have focused on a regime where the num- 
ber of excitations is much smaller than the total number 
of atoms, the laser driving is much stronger than the in- 
teraction between the Rydberg states and the atoms are 
distributed in a regular lattice. In this section we will 
test the validity of the first two approximations made in 
order to diagonalize the Hamiltonian, i.e. the Holstein- 
Primakoff approximation and the assumption that the 
dynamics is constrained to subspaces with equal excita- 
tion number. In addition, we will check the robustness 
of our results with respect to uncertainties in the atomic 
positions. This will bring us closer to the experimen- 
tal realistic situation with quantum uncertainty and/or 
finite temperature. 



A. Validity of Holstein-Primakoff and constrained 
dynamics approximations 

Throughout Section |III| where we performed the diag- 
onalization of the Hamiltonian (JlJ two approximations 
were made. First, the Holstein-Primakoff approximation 
TVb <C N allowed us to neglect all terms of the Hamil- 
tonian of higher than quadratic order. Second, we con- 
sidered that in the regime of strong laser coupling, i.e. 

^> V nn , the dynamics of the entire system is approx- 
imately driven by the Hamiltonian (J5|, which conserves 
the number of excitations. In order to test the validity 
of these approximations, we consider now the complete 
Hamiltonian (JlJ and rewrite it as 

H S» Eq + Hq + -ffpert 

where Eq is the unperturbed energy of the ground state 
Q and H is the unperturbed Hamiltonian ^ . The per- 
turbation to this Hamiltonian can be divided into three 
terms 



H, 



pert 



Hi + H 2 + Hup , 



where Hi and H 2 are given in ([6| and j7b, respectively 
and the third contribution Hup contains the terms of 



third and fourth order in a k that were neglected by the 
Holstein-Primakoff approximation 



JV 



H 



HP 



\ E { (44«fc + a W a k 



k=l 



2 E Vkm { a k a k a ln a k + a\a k a k a m 



k^m 

t t t t 

a k a k a k a m + a k a' m a k a k 



Let us now calculate the second order energy shifts 
caused by the three terms of H pcrt to the unperturbed en- 
ergies of the ground and first excited states of the Hamil- 
tonian Hq. These contributions vanish only in the limits 
|A| /fl -> 0, V nn /fl -> and N h /N ->■ 0. Here we con- 
sider that A = 0, and we calculate the contributions for 
a finite value of the other two ratios. 

We focus first on the corrections to the energy of the 
ground state. H pcrt only contains terms that couple 
states that differ by AN^ = 0, ±1, ±2. As a consequence, 
only the states of the first and second excited manifold 
([9]) and ( 10 ) contribute to the second order shift to the 



energy of the ground state, given by 



E™ = 



N 

E 



\(U\Hi\0)f 

Eq — El, 



E 



Eq — E<>.. 



Note that, since Hppp |0) = only Hi and H 2 contribute 
to the correction of the energy Eq 
yields 



(2) 



This correction 



E, 



(2) 



JV 

E 

i=i 



1 J^lAr l + 4|£ fe M fci 



40 + A 



(19) 



and depends only on the values of fi, N and the eigen- 
vectors and eigenvalues of the interaction matrix V 
given by M k i and A, respectively. The relative error 
Eq /Eq x 100 introduced by these shifts for a square 



lattice of side L = V N and several values of the ratio 
Q/V nn are given in Table [[] One can observe that the de- 
pendence with Q/V nn can be very well approximated to 
be proportional to (Q/V nn )~ , as it was already consid- 
ered in Section [Hi] Note that for n/V nn = 20 the shift 
of the energy of the ground state is below 0.6% for all 
sizes of the lattice considered, so that we can consider 
the unperturbed energy a very good approximation to 
the exact solution. The dependence with the size of the 
lattice has a more complicated form, although one can 
observe clearly that the perturbation becomes stronger 
for larger sizes N. This is expected since the number of 
states that are coupled to the ground state also increases, 
and so do the number of terms of the sum (19). 



We follow the same procedure to calculate the second 
order energy shifts of the energies of the manifold with 
one bosonic excitation. In this case, we also have to con- 
sider the coupling of the |1 Q ) states to the triply-excited 
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L 

ll/V 


3 


4 


5 


6 


7 


8 


9 


5 


5.0 


6.5 


7.5 


8.2 


8.7 


9.1 


9.4 


10 


1.2 


1.6 


1.8 


2.0 


2.1 


2.2 


2.3 


20 


0.31 


0.40 


0.45 


0.49 


0.52 


0.54 


0.56 


50 


0.049 


0.063 


0.072 


0.079 


0.083 


0.087 


0.090 



TABLE I: Error (in %) introduced by the second order energy 
shifts to the ground state energy Eq given in Q for A = 0, 
different values of Q/V nn and sides of the square lattice given 
by L = VN. 



L 


3 


4 


5 


6 


7 


8 


9 


5 


7.2 


8.1 


8.7 


9.1 


9.4 


9.7 


9.9 


10 


1.6 


1.8 


2.0 


2.1 


2.2 


2.3 


2.3 


20 


0.37 


0.44 


0.49 


0.52 


0.54 


0.56 


0.58 


50 


0.057 


0.069 


0.077 


0.082 


0.086 


0.089 


0.092 



TABLE II: Error (in %) introduced by the second order en- 
ergy shifts to the singly-excited state | ljv) energy E\ N given 
in Q for A = 0, different values of Q/V nn and sides of the 
square lattice given by L = yN. 



ones by H 2 and Hup: 



(2) 



y, |(0|gi|l a )| 2 >y | (2ij | Hi + H RP \l a ) 
^ E x -E Q ^ 



\~* W'i.i 

i<j<fe 



i<j 

\(3ijk\ H 2 + Hhp |1 q 



E\ a — E 2 . 
1 2 



Ei a - E 3ijh 



The final expression of the energy shift is in this case 
much more unwieldy and thus we will not include it here. 
As an example, we choose again the state of the single 
excitation manifold with highest energy, i.e. |1jv), that 
corresponds to the largest intensity of the transition from 
the ground state (the results are almost identical when 
choosing other states). The percentages of variation of 
the energy of this state due to the second order shift are 
shown in Table Il]for several values of f2/V nn and L and 
A = 0. Note that the values are similar to the ones of 
the ground state. The dependence on the laser strength 
follows again the expected (V na /Cl)~ behavior. The cor- 
rection to the energy due to the term 7?hp decreases with 
increasing number of sites N. This is expected, since the 
contribution of the terms of higher than quadratic order 
vanishes when N/Nb — > 0. On the other hand, the en- 
ergy shifts due to the other two terms of the perturbation 
Hamiltonian Hi and H 2 become larger when increasing 
the number of states (the size of the system) . The inter- 
play of these two opposing tendencies makes the depen- 
dence of the quantity E\' /Ex N x 100 less pronounced 
than in the case of the ground state. 



B. Uncertainty in the atomic positions, a/a 7^ 

We have considered so far that the atoms populate the 
harmonic oscillator ground state and that the atoms were 
infinitely localized, so that the width of the external wave 
function was negligible compared to the interparticle sep- 
aration, i.e. a /a — > 0. Let us now consider a small but 
finite value of this ratio given by the inevitable uncer- 
tainty of the atomic position due to the finite strength of 
the confinement and finite temperature of the system. 

We are particularly interested in how the intensities of 
the transitions from the ground state to the single excita- 
tion manifold I\ given by Eq. ( 14 ) change when consid- 



ering this uncertainty. To do so, we distribute randomly 
the position of the atoms around the center of each site 
of the lattice with a finite standard deviation a 7^ 0. As a 
consequence, the symmetry of the system is lifted, the se- 
lection rules discussed in Section |IV A| are not applicable 
anymore and thus more states of the single excitation 
manifold become accessible from the ground state. In 
Fig. [8] we show the resulting absorption profile after av- 
eraging over 10000 realizations up to a maximum value 
of a /a = 0.05 for a square lattice of 7 x 7 sites. We 
observe that as we increase the uncertainty in the posi- 
tion of the atoms the initial sharp lines (indicated by the 
red vertical lines) become broader and blue-shifted. The 
asymmetry of the shift can be explained as a consequence 
of the geometry of the lattice. Let us first consider the 
configuration with a = and that only one atom in the 
lattice can move from the center of the site. This atom 
can move in 6 possible directions, i.e. ±x, ±y or ±z. 
If it moves in the icy-plane, it will end up closer to at 
least one of its neighboring atoms. Since the interaction 
grows as 1/r 6 , the energy in this case will be increased. 
Only the configurations in which the atom moves in the 
±z directions will give rise to a lower interaction energy. 
In total 4 out of the 6 possibilities give rise to a larger 
energy, which explains the overall shift to the blue in the 
profile. When the atom is located in one of the bound- 
ary sites, the number of configurations that reduce the 
energy is 3 out of 6. As the number of sites on the lattice 
TV increases, the relative number of atoms on the bound- 
aries decreases and, as a consequence, the shift to the 
blue is more pronounced. This phenomenon is consis- 
tent with the observations obtained from the numerical 
simulations. 



VII. CONCLUSIONS AND OUTLOOK 

In this work, we have studied the excitation properties 
of a laser-driven gas of highly excited atoms confined to 
a two-dimensional lattice. In order to solve the system's 
Hamiltonian we have made two approximations. First, 
we have restricted our study to the states with a small 
number of excitations compared with the total number 
of atoms in the system. Moreover, we have focused on 
the regime where the strength of the laser is much larger 
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FIG. 8: a: Normalized absorption profile for the transition 
from the ground to the first excited states when the positions 
of the atoms in the lattice are randomly distributed around 
the center of each site with standard deviation a (in units of 
a). The red vertical lines show the transition lines for a — 0. 
The horizontal black line indicates the cut shown in b: for 
a — 0.025. c: The same cut for a — 0.05. The results are 
averaged over 10000 realizations. 

than the interaction between the atoms in the Rydberg 
state. After these approximations the resulting Hamil- 
tonian could be solved by diagonalizing a small matrix 
which allowed us to obtain the many-body eigenstates of 



the system and their corresponding energies. 

Moreover, we have described a protocol of how to ex- 
perimentally access the many-body states in the low- 
occupation sector of the spectrum by the use of an oscil- 
lating detuning. Here we have found that only the states 
with the same symmetry properties as the ground state 
of the system can actually be excited. We have shown 
how to map these collective many-body atomic states 
into non-classical single photon sources and observed that 
the interplay between the interatomic distance and the 
wavelength of the emitted photon changes dramatically 
the angular distribution of the photon state. In partic- 
ular, we have identified a regime in which the photon 
is emitted in a superposition of well-separated beams, 
whose number and relative direction depends strongly 
on the geometry (square or triangular) of the underlying 
lattice. This feature has potential applications in quan- 
tum information processing, especially to form quantum 
networks. 

Finally, we have considered a small uncertainty in the 
positions of the atoms on each site of the lattice due to 
the finite temperature and the quantum nature of the ex- 
ternal atomic wave function, which brings us closer to the 
actual experimental situation. The intensity of the tran- 
sition from the ground state to the first excited manifold 
has been calculated for varying values of the uncertainty 
in the position and we have found a systematic broaden- 
ing and blue-shift in the profile. 
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